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^~H A PNJL model is built, in which the Polyakov-loop potential is explicitly Zat^ -symmetric in order 

^^ to mimic a Yang-Mills theory with gauge group SU(A^c)- The physically expected large- A^c and 

Cn large-T behaviours of the thermodynamic observables computed from the Polyakov-loop potential 

fj are used to constrain its free parameters. The effective potential is eventually U(l)-symmetric when 

5^ Nc is infinite. Light quark flavours are added by using a Nambu-Jona-Lasinio (NJL) model coupled 

^^ to the Polyakov loop (the PNJL model), and the different phases of the resulting PNJL model 

are discussed in 't Hooft's large-A^c limit. Three phases are found, in agreement with previous 

CTN studies resorting to effective approaches of QCD. When the temperature T is larger than some 

Cn deconfinement temperature Td, the system is in a deconfined, chirally symmetric, phase for any 

^_^ quark chemical potential /i. When T < Td however, the system is in a confined phase in which 

fH chiral symmetry is either broken or not. The critical line T^{)j,), signalling the restoration of chiral 

O , symmetry, has the same qualitative features than what can be obtained within a standard Nc = i 

Oh 
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X 



PNJL model. 
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I. INTRODUCTION 

The structure of the QCD phase diagram is intimately related to our understanding of fundamental features of 
QCD, like for example confinement dynamics and chiral symmetry breaking, and to their interplay with in-medium 
effects like a nonzero temperature or quark density. This is the reason why a lot of effort is devoted to study this 
field, either on the theoretical side, to which the present work belongs, or on the experimental side through heavy- 
ion-collision experiments. Among the various effective frameworks used to study the QCD phase diagram (see e.g. 
the review [J), we will mostly focus on two of them: Polyakov-loop effective models for the pure gauge part of QCD, 
and the Nambu-Jona-Lasinio (NJL) model for the quark part. 

The Polyakov loop is defined as 

in which P is the path-ordering, g the strong coupling constant, Aq — A^Ta the temporal component of the Yang- 
Mills field, Ta the generators of the gauge algebra, and T the temperature. The integral runs on the compactified 
timelike dimension. The Polyakov loop is a precious tool to study the phase structure of a given Yang-Mills theory 
since {L{T,x)) = (7^ 0) when the theory is in a (de)confined phase 2J. Moreover, gauge transformations belonging 
to the center of the gauge algebra only cause L{T, x) to be multiplied by an overall factor. That is why it has been 
conjectured that the confinement/deconfinement phase transition in a Yang-Mills theory with gauge algebra q might 
be related to the spontaneous breaking of a global symmetry related to the center of g [3] ■ In the particular case of 
SU(A'c), deconfinement might thus be driven by the breaking of a global Z^r^ symmetry. The order parameter of the 
deconfinement phase transition should then be the traced Polyakov loop 

= ^Tr,L, (2) 

where the trace Tr^ is taken over the colour indices. The thermodynamic properties of pure gauge SU(3) QCD can 
then be studied by resorting to an effective scalar field theory where the potential energy density is Zs-symmetric, 
with e.g. the form ^ 



C/ = T*A 



-^l0P + ^H^ + ^( 



(3) 



The real coefficients hi can be fitted on lattice data. Various applications of this formalism can be found for example 
in [S]. Note that, in the following, and L will generally be indifferently called Polyakov loop. 
The NJL model is based on the Lagrangian [5] 

G 

Cnjl = q[^l^^^ - mg)q + - [{qqf + {qijsrq)^] , (4) 

where q is the quark field, niq the mass matrix, and t the Pauli matrices when an SU(2) fiavour symmetry is 
considered. The interaction terms are such that the Lagrangian is chirally symmetric. The NJL model is designed 
to model chiral symmetry breaking and study many related phenomenological problems; the interested reader may 
consult the review [7j for more information. In the original NJL model, fermions are not coupled to the gauge field: As 
shown in |5] , the coupling of this model to the Polyakov loop can be achieved by minimally coupling the quark field to 
a gauge field of the form A^^ = Aq S^q, that formally appears as an imaginary quark chemical potential. The so-called 
PNJL model resulting in this coupling has motivated a lot of studies devoted to the QCD phase diagram [9l [TO] . 
including cases with a nonzero magnetic field [IT] or nonlocal extensions [T^l [13] . 

The phase structure of the PNJL model at arbitrary Nc has been discussed in [Tl], as well as its large- A^c limit. One 
of the ingredients of this last work is to set A ex (A^^ — 1) in © so that the gluon potential has the correct scaling in 
Ac. The phase diagram that has been found at large- A'c is given in Fig.[T] The chirally symmetric but confined phase 
that appears for quark chemical potentials larger than about a third of the nucleon mass, fi > A/jv/3, can presumably 
be identified with the quarkyonic phase, that has been first proposed in 15 and further studied in [16 in particular. 

In the present work, we propose to re-build a PNJL model valid at large- A^, but in which the Polyakov-loop potential 
is explicitly Z^^ symmetric. A possible way to answer to this requirement is to introduce a term in (j)^" +(p*'^'= instead 
of the standard Z3-symmetric term (j)^ + cj}*^. Such a potential is proposed in Sec. |ll] and the corresponding PNJL 
model is written in Sec. |III[ Then, the issue of deconfinement and chiral symmetry restoration when varying T and fi 
are discussed in Sec. |IV| in 't Hooft's large- A'c limit. The obtained phase diagram and concluding comments are given 
in Sec. |Vj 
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FIG. 1. Phase diagram obtained by taking the large- A'^c limit {Nc — > (xj) of the PNJL model used in |14| . The solid lines signal 
first-order phase transitions. 



II. PURE GAUGE SECTOR 



A. Explicit Z]v^ -symmetry 

The simplest effective potential energy density depending on cj), defined in ([2]), and being explicitly Z^r^ -invariant 
has been proposed in [17] and reads 



*Na 



VgiT, N,, ^, 0*) - AiT, N,) \cj)\' + B{T, N,) \c^\' + C{T, N,) {(b''^ + <!>*''') 



(5) 



It is fornrally valid for any value of 0, but one may restrict oneself to |0| G [0,1] in a mean-field approximation. 
The above expression contains the basic blocks that could be expected to build a nontrivial theory: A mass term 
(|(/)p), an interaction term (|(/'|'*), and the term in (f>^^ + ip*^" accounting for the explicit ZTv^-symmetry THl. Terms 
scaling like \(j>\^, \(f>\^, ■ ■ ■ , l^]^'^^^, etc. could be added, but then the number of arbitrary functions would become 
too large to be efficiently constrained. Moreover, such higher-order terms would mostly be interaction terms that 
are already present in their simplest form in the |(/)|^ term. The expression (l5| is thus particularly convenient since 
it contains the minimal number of terms needed to perform the present study. The real coefficients A^ B, and C 
appearing in (l5| are functions of T and N^ , and their explicit form will be specified in the following. Note that 
0, which depends on T, N^, and x a priori, is here assumed to be independent of x. Beyond the polynomial form 
([5]), logarithmic shapes can actually be shown to emerge from a Haar integration on the gauge group in a strong 
coupling expansion. One can find such a form in [S], or for example in [19j, where a potential schematically given by 
U/T"^ = A{T)\(I)\'^ + B{T) In [l - Q\(j)\'^ + 4(03 _^ 0*3) _ 3|0|4j jg ^g^^ j^j. jv^ ^ 3 computations. Instead of computing 
a similar potential at arbitrary Nc, we keep the ansatz ([5| in the following; it is indeed particularly convenient for the 
calculations that are to be performed and still contains the Z^r^ symmetry we want to take into account. 

Various parametrizations of Za-symmetric potentials, fitted on pure gauge lattice data, have been proposed so 
far 00 HO]- Here, we are rather interested in obtaining an effective potential valid at large TVc, i.e. Nc> ^ aX least. 
As shown below, all these values of Nc will have in common that, in our approach, the "asymptotic" behaviour of Vg 
(i.e. values of |0| larger that the physical one minimizing the potential) will be driven by the Z^v^ -symmetric term. 
The present formalism will thus not be valid for Nc = 3 in particular, where this asymptotic behaviour is driven by 
the interaction term. The following expected qualitative behaviours have to be imposed in order to constrain the 
shape of the functions A, B, and C: 

• The pressure Pg = — min0(Vg) is proportional to NcT^ at large Nc and T in order to recover asymptotically 
the Stefan-Boltzmann limit of a free giuon gas. 



The norm, |0o|j of the optimal value of the Polyakov loop. 
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is A^c-independent at the dominant 



order, see the definition ([2|. The first corrections, scafing as l/Nc, are neglected in the present approach - more 



results on large- A^c features of Wilson and Polyakov loops can be found for example in [21]. |0o| 
confined phase, and > in the deconfined phase. Also, |0ol tends toward unity at very large T. 



in the 



There exists a critical temperature Tfj signalling a first-order phase transition, i.e. the potential must have two 
different minima whose depth changes with the temperature in order to modify discontinuously the localisation 



of the absolute minimum. At the critical temperature, |(/)o| = and 1/2 are two degenerate minima of Vg. This 
last value is chosen so that it will ensure a good compatibility between our model and existing lattice data but 
it has only to be nonzero in order to lead to a deconfined phase. T^ has to be seen as a typical value for the 
deconfinement temperature in SU(iVc) Yang-Mills theory since the deconfinement temperature appears to be 
A^c-independent up to corrections in 1/N^ 22, 23j. 

Obviously, the above constraint does not apply to Nc = 2, where the transition is of second-order. This is not 
problematic since we eventually look for a model valid at large- A^c- Moreover, the value \(j)o\ = 1/2 may not be the exact 
value of the Polyakov loop in T^: Recent lattice results find it to be around 0.4 f24j, while a more recent renormalization- 
group-based approach leads to values closer to 0.6 for the Polyakov loop at the deconfinement temperature ;23j. The 
value 1/2 then appears to be relevant because it falls in the typical range of the existing results and because it 
simplifies the calculations performed in the following. 

The above constraints are actually satisfied by the following Lagrangian 



Vg = N^T^ a{T) 
where 



101' - m' + ^-^^^mrf - 1](0^= + 0*^^) 



(6) 



a(r)>0, liT)>^, l{Td)^l, dTl{T)>0, ?(oo) = l. (7) 

Explicit forms of a{T) and l(T) will be given in the next section. All these conditions are required in order to have 
the existence of 2 degenerate minima and the correct behaviour of the Polyakov loop in the mean field approximation. 
The potential ^ has the following absolute minimum: (j)o{T < Trf) = and (f)o{T > Td) = |</)o(T)| e^'^'^'/A'c^ where 
k = 0, . . . ,Nc— 1 and where |0o(2^)| is a solution of 

1 - 8\MT)\' + KTf-"^ [SliTf - 1] |0o(r)r-2 = 0. (8) 

It is straightforwardly checked that 

\MT)\^l{T) (9) 

actually solves (Is]). 

A more compact expression for the optimal value of the Polyakov loop is thus 

(t>o = ;(r) 62^'^'=/^'= e(T - r^), (lo) 

where Q is the Heaviside function. As seen from (Is]), \(j}o\ only depends on T as required. 

Restricting ourselves to the values = |0| e^*^'^/"^^ we get at the limit A'c — > oo a quite simple shape for the effective 
potential ([6]), namely 

^' -'^^ = a(r)|0|2(l-4|0|2) \<f>\<l{T), 



^-Hoo \(t)\>l{T). (11) 

Hence, a U(l) invariance is recovered at infinite Nc as a limiting case of the Zjv^ -symmetry. The schematic evolution 
of the large- A"c limit of Vg with the temperature is plotted in Fig. [21 the behaviour (111 is readily observed, as well as 



the change of global minimum vclT — T^. Finally, the large- A'j, limit of the pressure reads 

Pg{T,Nc) = NlT^a{T)l{Tf [Al{Tf - l] . (12) 

Provided that l{oo) — 1 according to the large-T behaviour of the Polyakov loop, Pg would tend toward the Stefan- 
Boltzmann limit for a free gluon gas if a(oo) = 7r^/135. 

B. Numerical data 

The function 1{T) is constrained by the relations ([7|) in order for the structure of the potential and its evolution 
with the temperature to have the required behaviour. Moreover, 1{T) is equal to the norm of the Polyakov loop as 
soon as T > Td- Those physical constraints are not sufficient to write down an explicit expression for 1{T). A possible 
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FIG. 2. Schematic evolution of the effective potential ( 11 1 versus the temperature (solid lines) 
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FIG. 3. Norm of the Polyakov loop minimizing the potential m versus the temperature in units of T^ (solid line). The function 
1{T) (dashed line) and the norm of the Polyakov loop computed in pure gauge SU(3) lattice QCD (points) have been added for 
comparison. Lattice data are given for temperatures lower than 2.4 Td; data are taken from [24j . 
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way of proceeding, that we choose here, is to fit 1{T) on available lattice computations of the Polyakov loop in pure 
Yang-Mills theory. To our knowledge, large-iVc values have not been obtained so far, but accurate SU(3) ones have 
been computed in |21]. Since the Polyakov loop should not depend on Nc at the dominant order, it is relevant to fit 
1{T) on SU(3) data; the ad hoc form 



1{T) =0.74-0.26 tanh 



2.10 



0.60 



T 



(13) 



leads to a satisfactory parametrization of the results of [25 as it can be seen in Fig. [SJ It is also worth noting that 
Fig. [2]has been obtained using the form (13) for 1{T). 

It is important to remark at this stage that the calculations we will perform are done in the mean-field approximation. 
In this scheme, the Polyakov loop is always lower than 1: values larger than 1 are due to quantum fluctuations and are 
de facto beyond the mean- field treatment. That is why we have restricted our fit to lattice data lower than unity (T < 
2.4 Td). We miss the overshoot due to quantum fluctuations, but we stay coherent with the mean-field approximation, 
and reach moreover Z(oo) = 1. As a consequence, our results should be mostly trusted below T < 2.4 Td but this is 
not a fiaw since, in the following, we will be concerned with the phase structure of the theory and no phase transition 
will appear at energy scales above this upper limit. 

The positive-definite function a{T) is only present as an overall factor in Vg, so it does not come into play in the 
qualitative features of the effective potential. However, it is relevant in view of reproducing the absolute value of the 
pressure in pure gauge QCD, for which lattice data are known at Nc = 3, 4, 5, 6, 8 and cx) through an extrapolation 
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FIG. 4. Large- Nc pure gauge pressure computed from Eq. (121 and normalized to N^T"^ (solid line). The corresponding lattice 
data, taken from I 



are plotted for comparison in the case A'^c = 3 (gray points) and Nc -^ oo (black points). 
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The empirical choice 



a{T) = 



1{TY 



n 
135 



0.029 
\n{T/Td + 1.5) 



(14) 



leads to a good agreement between the lattice data of [25, and formula ( [l2| , as shown in Fig. [4] Notice that the value 
a(oo) is such that the Stefan-Boltzmann limit is reached at large temperatures. 

It is worth summarizing what has been done at this stage. Starting from Lagrangian ([5]), we have shown that 
the three arbitrary functions of T and Nc it contains can be strongly constrained by demanding that the averaged 
Polyakov loop and the pure gauge pressure have a relevant behaviour in Td, at large T, and in the large- A^c limit. 
Explicit forms for the two remaining unconstrained functions of T can then be found by asking the present model to be 
in agreement with current pure gauge lattice data. One is finally left with a fully determined Lagrangian with explicit 
ZjVc symmetry at finite N^ and U(l) symmetry in the large-iVc limit. This Lagrangian is obviously not predictive 
concerning the thermodynamics of the pure gauge sector, just as previously used Lagrangians like ([3]). However, its 
knowledge is a necessary step in view of making predictions concerning the quark sector, whose inclusion is discussed 
in the next section. 



III. PNJL MODEL 



As shown in [8] , a minimal coupling of the N JL Lagrangian ( 4]) to a gauge field of the form A^ = Aq (5^o makes 
eventually appear the Polyakov loop in the quark grand potentiaL In the mean field approximation, one is led indeed 
to the quark potential [8 



N,Nj 2g J (27r)3 

T 
Ep + — Tr^ln 



(15) 



l + Le-(^--'^'/^ 



N, 



-Tr^ln 



l + Lte-(Sp+P)/T 



where the Polyakov loop L has been defined in (fTl). In the above equality. 



Ep = ^p^ + (m, - a)2 
is the quark dispersion relation, with niq the quark bare mass and a related to the chiral condensate as follows 

a = G{qq). 



(16) 



(17) 



The coupling G has to scale as {Nc Nf) ^ in order for the potential ( 15 ) to scale as Nc Nf, so it is convenient to define 
the coupling g as 



g = GNcNf. 



(18) 



Although Nf is a priori arbitrary, our results are mostly valid at Nf = 2. For higher values of Nf indeed, the axial 
anomaly (not present in this formalism) should be taken into account in order to get a reliable model. In what follows, 
Nf = 2 will be implicitly understood, although we keep the notation Nf so that the quark contributions appear more 
clearly. 

Since the pure gauge part of the potential only involves the traced Polyakov loop (j), it is interesting to express Vq in 
terms of cj) rather than L. Terms of the form Trcln [1 + z L] can be expressed as functions of TicL oc 0, TtcL^, TrcL'^, 
. . . through a Taylor expansion. A possible way of proceeding is to expand the quark potential at the first order in 
L. This eventually leads to formulas in which only (f) appears in Vq |14j . This scheme has the advantage of being 
independent of the parametrization of L. Here we adopt an inequivalent procedure. As a first step, we notice that 
there exists in general a gauge in which the Polyakov loop L is a diagonal element of SU{Nc): 



X = diag(e*^\ e'' 



pltfNc-l 



ENc — 1 a 



(19) 



The Nc — I parameters 6j are real so that L^'^L = 1 and deti = 1 as demanded for an SU(A^c) element. In the special 
case of Nf, = 3, there is a one-to-one correspondence between the parameters 6i, O2 and the Polyakov loop degrees of 
freedom (/), (jf . This is not the case at large N^ however, where the number of independent parameters in the Polyakov 
loop goes to infinity. As a consequence, an exact computation of the color traces appearing in ( 15 ) is not possible 



unless simplifying assumptions are made. As a second step to reach this goal, we propose the following ansatz: 



L = diag(e*^...,e*^l,e-*' 



diag(e'^...,e*^e- 
^^ ■' '■~- 



(Ar,-l)/2 



odd - Nr 



even — N^- 



(20) 



iVc/2 



Na/2 



It reduces to the mean-field parametrization of [S] at TVc = 3, but the price to pay is that the number of degrees of 
freedom in L is drastically reduced to a single real parameter 9. It is then readily computed that 



= ^[l + (iV,-l)cos( 



odd - Nc 
even — Nr 



(21) 



by using of the ansatz (201 in ^. We thus have an asatz that "looks like" the SU(3) case and that reduces to (f> — cos 6 
at large- A'c. 

Moreover, one can compute that 



Tr^ln 



l + Le-(^--^)/^ 



= Indetc 
iV, -1 



l + Le-(^--^)/^ 
In 



Nr-l 



-In 



1+e 



-(Bp-M)/T 



odd - Nr 



(22) 



Nr 



In 



1 + 2,^e-(^''-^)/^ + c-2(Bp-M)/T 
even — Nc, 



and, taking into account a cutoff for the momentum integration of the vacuum term, one finally arrives at the quark 
potential, whose large-A^c limit is given by 



ujq{fi,T,a,( 



Vq{fi,T,a,cl)) 

NcNf 

2 1 /-A rp 

27r2 



a- 1 



dpp Ep 



dpp" {in [1 + 20e-(^--^)/^ + e-2(Sp-M)/2 



(23) 



This last potential formally reduces to the genuine NJL potential once </> = 1, as observed in previous studies p.fH]. 
The total potential of the large- A^c PNJL model under study is finally given by 



V(/i, T, a, 0) = NlLOg{T, 0) + Nc Nf a;,(/x, T, ct, , 



(24) 



In the confined phase, where (/) = 0, one observes a term in In [l -I- e 



-2{Ep-^,)/T^ 



in the potential (23), so it could 



be tempting to associate such a term with diquark degrees of freedom in the confined phase. However, in the limit 



T ^- 0, one exactly recovers the zero temperature NJL potential, expressed in terms of the quark degrees of freedom. 
So the confined degrees of freedom are still quarks in the present approach. Similarly, baryonic degrees of freedom 
(states with Nc quarks) are not included in the present formalism; this would require further extensions of the present 
approach that are beyond the scope of the present study. 

IV. PHASE DIAGRAM AT LARGE A^,, 

In 't Hooft's large- A'c limit, the number of quark flavours stays finite and V is dominated by the gluonic contribution. 
Consequently, when Nc becomes infinite, the optimal value 0o can be found by minimizing ujg only. According to 



(10 1, the large- A'c solution (i.e when N^ is infinite) reads 

MT)^l{T)e{T-Td). (25) 

The physical value of a, denoted (Jq and depending on T and /i, is then such that it minimizes iOq{T, ii,a,(f)Q{T)). 
Ldg does not depend on a. Since a oc (qq), chiral symmetry is present when ao — and broken when ctq ^0. As a 
consequence of (25), the deconfined phase appears as soon as T > T^, independently of the value of fi: As pointed 
out in [15' , quarks have no influence on the deconfinemcnt phase transition at large- A'c because of the suppression of 
internal quark loops in this limit. 

As a consequence of the large- A^c limit, the confined/deconfined phases are straightforwardly identified in our model. 
The situation is less simple as far as chiral symmetry is concerned; numerical computations are needed. As a first 
step, the parameters of the model have to be fixed. The values 

niq = 5.5 MeV, g = 60.48 GeV~^ A = 651 MeV, Ta = 270 MeV, (26) 

used in the PNJL study [9], will be taken in the following also. The first three parameters have been fitted so that 
the zero-temperature pion mass and decay constant are reproduced within the standard NJL model with iVc = 3 and 
Nf = 2 [71 [55]. Td is a typical value for the deconfinement temperature in SU(A'c) Yang-Mills theory. 



Using the parameters (26 1, the optimal value ctq can now be computed for any couple {fj,,T), and can be linked to 



the quark condensate thanks to (17) 



{qq){^l,T)^^^^ao{^i,T). (27) 

9 

In the limit where T and /i both tend toward zero, we get 

lim (qq) {fi, T) = -N^ Nf 5.29 10^ MeV^ (28) 

corresponding to a quite common value of —(317 MeV)'^ for TVc = 3 and Nf = 2. 

The large- A'c chiral condensate versus the temperature is plotted in Fig. [5] for some values of the quark chemical 
potential. The most salient feature of this plot is the simultaneity of the first-order deconfinement phase transition 
and of the restoration of chiral symmetry through a first-order phase transition occurring at T^ = T^. However, 
when n/Td > 0.8 (/i > 200 MeV), the quick decrease of the chiral condensate suggests a progressive restoration of 
chiral symmetry through a crossover at temperatures smaller than T^. As shown in [5], the crossover temperature 
can be computed thanks to the determination of the peak position in the dimensionless quark susceptibility reading, 
at large- A'c, 

o- q\(j—ao 

We have chosen to follow the definition of [8^ for the quark susceptibility, although a more standard way of defining 
the susceptibility is rather d^ ujq{T,pL), see e.g. \\Z\. In both cases, a peak in the quark susceptibility signals a phase 
transition . 

A plot of Xqq{T,fJ-) for some values of n/T^ is given in Fig. p] Several observations can be made by observing 
this figure together with Fig. [5] First, the peak of the quark susceptibility is located in Td when /i/T^ < 0.79; 
this corresponds to a first-order-type chiral symmetry restoration in the deconfined phase. The point (0.79, 1) x Td 
actually corresponds to a triple point in the (/i, r)-plane: At large /i the peak of Xqq is located below Td - a larger 
fi corresponds to a lower peak position -, leading to the existence of a confined phase in which chiral symmetry is 
progressively restored through a crossover. A careful look at ctq actually shows that the chiral phase transition below 
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FIG. 5. Chiral condensate at large- A^c {Nc — >■ cxa), normalized to its zero temperature value, versus T in units of Td, and plotted 
for n/Td = 0, 1, and 1.11 (solid lines). The optimal value of the Polyakov loop is also plotted (dotted line). 
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FIG. 6. Dimensionless quark susceptibility (29 1 versus the temperature in units of Td (solid lines) with, from left to right, 
fi/Td =1.24, 1.15, 1.04, 0.93, 0.79. XqqiTtO) is also plotted for completeness (dashed line). Computations were done for 

Nc — > oo. 



Td becomes of first order when fJ,/Td > 1.24: There exists thus a critical-end-point that we find to be (1.23, 0.26) x Td in 
the (/i, r)-plane. Apart from the susceptibility, the position of the chiral phase transition could have been alternatively 
determined by computing the zero of d^aoiT, /i). We have checked that the chiral temperatures computed using that 
method agree with those computed with those computed from the peak in the susceptibility up to 5%. For an 
exploratory study such as the present one, this agreement is satisfactory. 

Gathering all these observations, the phase diagram of our model in the (/i, r)-plane can be established; it is shown 
in Fig. It] The three phases we find correspond to those found in [M] , see Fig. [11 but the structure of the chiral phase 
transition is a bit more involved under T^: The chemical potential at which chiral symmetry is restored now depends 
on r, and there exists a critical-end-point at large enough fj,. Although the deconfining phase transition corresponds 
to what is expected in the large- A^c limit of QCD from generic arguments |15j . the critical line Ty^{^) we find under 
Td quite resembles to what can be observed within previously known Ac = 3 PNJL studies (H |9] . The similarity 
between our way of including the Polyakov loop in the NJL model and the way of [5| - our ansatz is a straightforward 
generalization of the one used in this last work - might actually be at the origin of the similarities between the phase 
diagrams we find. The same reason, combined to the fact that we chose for our parameters values fitted on the SU(3) 
case, might explain why the values we find for e.g. the chiral condensate are similar to those of [8]. 

We notice that, at large but finite values of N^ the full potential ( p4| has to be minimized and quark contributions 
(presumably in 1/Ac) will cause the Polyakov loop to be different from 1{T). Hence, the chiral condensate will also be 
modified, and the whole phase diagram will be affected. We nevertheless choose here to focus on the large- A"c limit 
of the model, since it has been designed to be relevant in this limit mostly. 
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FIG. 7. Phase diagram of the large- A^c PNJL model (24 1 with explicit Zat^ symmetry, obtained for Nc — )■ oo. The solid lines 
denote first-order phase transitions while the dashed line denotes a crossover. The triple point (0.212, 0.270) GeV and the 
critical end-point (0.335, 0.063) GeV have been also plotted. The end of the lower curve is reached at (0.343, 0) GeV. 



CONCLUSIONS 



Effective "Polyakov-loop-based" approaches have proven to be a relevant tool in view of modelling the thermody- 
namic properties of pure gauge QCD. The traced Polyakov loop is then the order parameter associated to confinement, 
itself seen as correlated with a global center symmetry, Z^^ when the gauge group is SU(A'^c)- Following a suggestion 
made in [T^, an explicitly Zat^ -symmetric potential involving the traced Polyakov loop has been built and leads to: 
A first-order phase transition at large- N^ a gluonic pressure scaling as N^, and an TVc-independent optimal value for 
the Polyakov loop. The coupling of the pure gauge sector to light quarks has then been performed within a PNJL 
approach. Thanks to a particular ansatz for the Polyakov loop, the quark potential is such that it only involves the 
traced Polyakov loop that appears in the pure gauge potential. It has to be said that the resulting PNJL model is 
designed to be relevant in the large- iVc limit only. It is moreover inequivalent to the model proposed in [T3j since the 
colour traces that have to be performed in the derivation of the PNJL potential are carried differently. 

The different assumptions that were made lead to a large- TVc phase diagram that compares favourably with pre- 
vious approaches. At any /i, it shows a deconfined, chirally symmetric, phase above the deconfinement temperature 
{Td =270 MeV). The deconfinement phase transition and the restoration of chiral symmetry are found to be simulta- 
neous first-order phase transitions, in agreement with |14[ 115] . At temperatures lower than Td, thus in the confined 
phase, there is a critical line T^{fJ-) separating a phase with broken chiral symmetry at small /i and a chirally sym- 
metric phase at large fi. The phase transition is found to be a crossover from the triple point (0.212, 0.270) GeV 
to the critical-end-point (0.335, 0.063) GeV. It is then of first order until the boundary (0.343, 0) GeV is reached, 
corresponding to the estimate fi sa Mn/3 [TS]. It is worth saying that a confined, chirally symmetric, phase has also 
been found by solving Schwinger-Dyson equations at nonzero /Lt in [27], and that evidences for the existence of such 
a phase has been found in Coulomb gauge QCD calculations [25]. The present model's prediction of three different 
phases, among which a chirally symmetric phase at baryonic chemical potentials around the nucleon mass thus seems 
a reliable one. The structure of the critical lines appears however to depend quite strongly on the way the Polyakov 
loop L is handled: The comparison of Figs, [l] and [7] shows how a different treatment of the Polyakov loop within 
similar PNJL-based approaches can affect the phase diagram. It is also worth noting that an accurate description of 
the confined but chirally symmetric - quarkyonic at large Nc - phase would require a more detailed approach that 
the simple PNJL model presented here. In particular, the explicit inclusion of baryonic degrees of freedom should be 
performed, but this is a task that we leave for subsequent studies. A first interesting attempt has for example been 
made in |29| . where a constituent approach is used to model the baryonic matter at large Nc- 

Recently, an effective model has been proposed [3D], in which the Lagrangian involves a linear sigma model for 
the quark part and a dilaton-like effective potential for the gluon part. Such a dilaton-like potential is designed to 
mimic the breaking of scale invariance through a nonzero value of the gluon condensate (the dilaton field) . The phase 
diagram which is found is nearly identical to the one of Fig. [7| excepted that the deconfinement phase transition is 
replaced by a restoration of scale invariance. Both seem then to be linked, as suggested in [30_. 

Finally, we remark that the large- TVc limit of the proposed pure gauge effective potential has a U(l) symmetry, which 
emerges as the limit of a Zjv^ symmetry. This large- A^c effective potential could be used in an approach where the 
traced Polyakov loop is allowed to depend on the position, typically via a Lagrangian of the type £ w dfj,(j)d'^4>* —Vg. 



11 



Of particular interest would then be to search for locaHsed, solitonic, solutions of C: One could then take advantage 
of the fact that finding solutions of a complex scalar field theory with a U(l) invariance is a topic that has attracted a 
considerable attention, mostly since Coleman's work on Q-balls [3T], and for which many results are already available. 
We hope to present such a study in future works. 
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